Towards an Asymptotic-Safety Scenario for Chiral Yukawa Systems 
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We search for asymptotic safety in a Yukawa system with a chiral U(JVl)l®U(1)r symmetry, 
serving as a toy model for the standard-model Higgs sector. Using the functional RG as a nonper- 
turbative tool, the leading-order derivative expansion exhibits admissible non-Gaufiian fixed-points 
for 1 < JVl < 57 which arise from a conformal threshold behavior induced by self-balanced boson- 
fermion fluctuations. If present in the full theory, the fixed-point would solve the triviality problem. 
Moreover, as one fixed point has only one relevant direction even with a reduced hierarchy problem, 
the Higgs mass as well as the top mass are a prediction of the theory in terms of the Higgs vacuum 
expectation value. In our toy model, the fixed point is destabilized at higher order due to massless 
Goldstone and fermion fluctuations, which are particular to our model and have no analogue in the 
standard model. 



I. INTRODUCTION 

The Higgs sector is a crucial building block of the stan- 
dard model of particle physics and parameterizes the 
masses of matter fields and weak gauge bosons. This 
successful parameterization goes along with two prob- 
lems of the standard model which give rise to the belief 
that the standard model should be embedded in a larger 
fundamental framework: the hierarchy problem and the 
triviality problem. These problems have initiated many 
further developments. 

The triviality problem renders the theory ill-defined 
from a fundamental point of view, since it inhibits an 
extension of the standard model to arbitrarily high mo- 
mentum scales P, S 0, 0, S @ ■ The scale of the max- 
imum ultra-violet (UV) extension Auv.max induced by 
triviality is expected to be related to the Landau pole of 
perturbation theory. Of course, the Landau pole, i.e., the 
divergence of a perturbative running coupling at a finite 
UV scale, in the first place signals the breakdown of per- 
turbation theory. Near the Landau pole, nonperturbative 
physics can set in and severely modify the picture. So, 
a study of the triviality problem and the existence of a 
finite Auv,max therefore requires a nonperturbative tool. 
This has also inspired a number of lattice investigations 
of scalar and chiral Yukawa systems @, B S 03 ■ 

Whereas the triviality problem is a true conceptual 
problem of the standard model, the hierarchy problem, 
i.e., the possible existence of a huge difference between 
the electroweak scale and an underlying fundamental 
scale such as the GUT or the Planck scale, is only a 
problem of unnaturally fine-tuned initial conditions (in 
this case of the mass parameter in the Higgs potential at 
the underlying scale). 

Solutions to these two problems are often sought by 
introducing field theories with new degrees of freedom 
or higher symmetries or by going beyond quantum field 
theory. The guiding principle of this work is more con- 
servative, as we intend to identify (or rule out) possi- 
ble solutions within quantum field theory essentially us- 
ing the same or very similar degrees of freedom of the 
standard model. Quantum field theory indeed offers 



a framework for such solutions in terms of Weinberg's 
asymptotic safety scenario [HI, [T3, EH, which has al- 
ready been investigated in a variety of models ranging 
from four-fermion models [lil . Il5l . Il6j |. simple Yukawa 
systems 17 1, nonlinear sigma models in d > 2 11811 . and 
extra-dimensional gauge theories to gravity [20L l2l| . 
For the asymptotic-safety scenario to apply, a fixed point 
of the renormalization group (RG) flow in the space of 
couplings has to exist. If the system sits on an RG tra- 
jectory that hits the fixed point in the UV, the UV cutoff 
can safely be taken to infinity, and the theory can remain 
valid to arbitrarily short distance scales. 

In fact, a suitable fixed point has recently been iden- 
tified in a Z2-symmetric Yukawa system [17| . serving as 
a simple toy model for the top-Higgs sector of the stan- 
dard model. Moreover, a generic mechanism inducing 
such a fixed point has been proposed which relies on 
a conformal behavior of the Higgs vacuum expectation 
value (vev). This conformal- vev mechanism is generated 
by a dynamical self-balancing of bosonic and fermionic 
degrees of freedom in the UV. In the simple Z2 model, 
this balancing actually occurs only for small fractional 
flavor number Nf < 0.3. This observation is a strong mo- 
tivation to consider more realistic Yukawa systems and 
to explore the potential of the conformal- vev mechanism 
for the UV problems of the standard model. Another 
scenario to circumvent the problem of triviality in such 
a Z2-symmetric Yukawa system in the same framework 
occurs by coupling the system to a gravitational back- 
ground [221 ] . 

This work is devoted to an investigation of a chiral 
^(A^l)l ® U(1)r Yukawa model, serving as a more so- 
phisticated toy-model for the Higgs sector of the stan- 
dard model (or of a GUT-like theory). The model is 
designed to have a left-handed chiral sector as is typical 
for the Higgs sector of the standard model. The num- 
ber of left-handed fermions N-^ is left as a free parameter 
in order to study the dependence of a potential fixed 
point on the varying numbers of degrees of freedom. As 
will be demonstrated by our analysis, this chiral struc- 
ture facilitates a more boson-dominated Higgs vacuum 
expectation value which is a prerequisite for conformal- 
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vcv mechanism to work. Using the functional RG as a 
nonperturbative approach, we systematically search for 
the existence of non-GauBian interacting fixed points of 
the RG flow which could allow for an extension of the 
model to arbitrarily high momentum scales and render 
the system asymptotically safe. 

A particularly attractive by-product of an asymptot- 
ically safe Higgs sector is given by the fact that the 
number of physical parameters is dictated by the num- 
ber of RG-relevant directions at the fixed point and is 
thus an inherent property of the model. As this num- 
ber can actually be smaller than the corresponding one 
at the Gaufiian fixed point (defining the "perturbativc" 
standard model), asymptotic safety can lead to a reduc- 
tion of physical parameters and hence have more pre- 
dictive power. For instance, the interacting fixed point 
in the Z2-invariant toy model has one relevant direction 
less than the perturbative fixed point. As an immedi- 
ate consequence, the value of the Higgs mass becomes 
a prediction once the Higgs vev and the top mass are 
fixed. For a fixed point of our chiral model discussed in 
this work, even the top mass can become a prediction, 
demonstrating the predictive and constraining power of 
the asymptotic safety scenario. 

A crucial question for all nonperturbative techniques is 
the systematic consistency and reliability of the results. 
In this work, we compute the RG flow of the model in a 
systematic derivative expansion of the effective action. 
This expansion is controlled if the momentum depen- 
dence of full effective vertices takes only little influence 
on the flow. A direct means for measuring this influence 
is the size of the anomalous dimensions rj of the fields, 
since next-to-leading order contributions couple to the 
leading-order derivative expansion only via terms ~ rj. 
Monitoring the size of rj thus gives us a direct estimate 
of the reliability of our results. Whereas the anomalous 
dimensions at the fixed point of the Z2 model were indeed 
found to be small, the anomalous dimensions at the fixed 
point of the present model can become large. Hence, the 
results within the present model have to be taken with 
a grain of salt. The reason for the difference between 
the two models lies in the existence of Goldstone bosons 
as well as massless fermions in the present model which 
contribute dominantly to the anomalous dimensions. As 
these massless modes are not present in the standard 
model, we expect that the leading-order derivative ex- 
pansion (where 77 = 0) of the present system can serve as 
a model for the Higgs- Yukawa sector of a more realistic 
gauged version. From this viewpoint, the essential idea 
to have a non-Gaufiian fixed point works and provides us 
with a highly predictive theory. 

This paper is organized as follows: in Sect, [ill we sum- 
marize the essence of the conformal-vev mechanism and 
introduce the concept of asymptotic safety. Sect. |TTT] dis- 
cusses the nonperturbative construction of the effective 
action in terms of the functional RG. The fact that the 
resulting flow equations do not support asymptotic safety 
in the symmetric regime is briefly elucidated in Sect. [TVl 



In Sect. El the regime of spontaneous symmetry breaking 
is explored to leading order, revealing non-Gaufiian fixed 
points with physically appealing properties. The result- 
ing predictive power of the asymptotic-safety scenario is 
described in Sect. I VII Sect . I VIll summarizes the problems 
of the present toy-model occurring at next-to-leading or- 
der in the derivative expansion. Conclusions are drawn 
in Sect. IVIIII Technical details on the derivation of the 
flow equations are given in the appendices. 

II. AN ASYMPTOTIC-SAFETY SCENARIO 
FOR YUKAWA SYSTEMS 

A. Conformal vacuum expectation value 

The model building of the present work is strongly 
motivated by qualitative considerations about the loop 
contributions to the running of the dimcnsionless version 
of the bosonic field expectation value v in the regime of 
spontaneous symmetry breaking (SSB) where v > 0. Our 
central idea is that the contributions with opposite sign 
from bosonic and fermionic fluctuations to the vev can 
be balanced such that the vev exhibits a conformal be- 
havior, v = (ip) ~ k. Here, k is the scale at which we 
consider the couplings of the system. The dimcnsionless 
squared vev k = ^v 2 /k 2 has a flow equation of the form 

v 2 d 
dtK. = dt-rriT = —2k + interaction terms, dt = k — . (1) 
2k z dk 

If the interaction terms are absent, the Gaufiian fixed 
point K = is the only conformal point, correspond- 
ing to a free massless theory. If the interaction terms 
are nonzero, e.g., if the couplings approach interacting 
fixed points by themselves, the sign of these terms de- 
cides about a possible conformal behavior. A positive 
contribution from the interaction terms gives rise to a 
fixed point at k > which can control the conformal 
running over many scales. If they are negative, no physi- 
cally acceptable conformal vev is possible. Since fermions 
and bosons contribute with opposite signs to the inter- 
action terms, the existence of a fixed point > cru- 
cially depends on the relative strength between bosonic 
and fermionic fluctuations. More specifically, the bosons 
have to win out over the fermions. 

In Fig. [T] (taken from [l7| ) , we sketch various options 
for the flow of the dimensionless squared vev k. The 
solid line depicts the free massless theory with a trivial 
Gaufiian fixed point at k = . If the fermions dominate, 
the interaction terms are negative and the fixed point is 
shifted to negative values (being irrelevant for physics), 
cf. dotted line. If the bosonic fluctuations dominate, the 
k flow develops a non-Gaufiian fixed point at positive 
values k* > 0. This can support a conformal behavior 
over many orders of magnitude, cf. dashed line. This 
fixed point is UV attractive, implying that the vev is a 
relevant operator near the fixed point. If the interaction 
terms are approximately k independent, the slope of dtK 
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FIG. 1: Sketch of the flow of the dimensionless squared Higgs 
vacuum expectation value k as described in the text below. 



near the fixed point is still close to —2, corresponding 
to a critical exponent O ~ 2 and a persistent hierarchy 
problem. An improvement of "naturalness" could arise 
from a suitable k dependence of the interaction terms 
that results in a flattening of the k flow near the fixed 
point, cf. dot-dashed line. Whether or not this happens 
in a specific model is a prediction of the theory which 
needs to be deduced from the theory by nonperturbative 
methods. 

In this work, we construct a model with standard- 
modcl-likc symmetries along this line of research. We 
introduce jVl left-handed fermion species ( a 
{1, ./Vl}) and one right-handed fermion ipn, as well 
as TVl complex bosons cf> a . All fields live in the funda- 
mental representation of the left-handed chiral symme- 
try group U(JVl)l- The Yukawa coupling is then realized 
by a term h('ipB,<j) a ' ip® — ^^""tpn) . This chiral Yukawa 
system mimics the coupling between the standard-model 
Higgs scalar and the left- and right-handed components 
of the top quark, also involving Yukawa couplings to the 
left-handed bottom (for N-^ = 2) and further bottom- 
like quarks (for > 2) in the same family. If the scalar 
field develops a vev upon symmetry breaking, the top 
quark acquires a Dirac mass, whereas the bottom- type 
quarks remain massless (similar to neutrinos in the stan- 
dard model). For = 2, we ignore the Yukawa coupling 
~ '0Re a &</' a ^L m our model, which provides for a mass 
term for the bottom quark in the standard model, since 
it does not generalize to other N^. 

In comparison to left-right symmetric models, this 
model has an interesting new feature concerning the rela- 
tive weight of the boson interaction terms contributing to 
the renormalization of the dimensionless vev k (sec fig- 
ure^: diagrammatically speaking, the inner structure 
of fermion/boson components in the fermion loop for a 
specific choice of external boson legs is fully determined. 
On the other hand, the boson loop obtains contributions 
from all Ni, components and so is linear in N-^. In this 
way, -/Vl serves as a control parameter for boson domi- 
nance and for the potential existence of a non-Gaufiian 
fixed point. Already at this qualitative level of the discus- 
sion, it is worthwhile to stress that the standard model 
has such a left-right asymmetric structure which can sup- 




FIG. 2: Loop contributions to the renormalization flow of 
the vev. The left loop involves only inner boson lines. The 
vertex A2 allows for a coupling between all available boson 
components. This implies a linear dependence on 7Vl for the 
renormalization of the boson contribution, see below. On the 
right panel, we depict the corresponding fermion loop contri- 
bution. The incoming boson 4> a fully determines the structure 
of the fermion loop and does not allow for other left-handed 
inner fermions than ip®, inhibiting an TVl dependence of the 
algebraic weight of this loop. 

port the conformal-vev fixed point. 

B. Flow equation and asymptotic safety 

In the asymptotic-safety scenario of a quantum field 
theory, the microscopic action to be quantized is a priori 
unknown. A construction of a renormalizable theory is 
possible if a suitable fixed point exists in the space of 
all possible action functionals, i.e., in theory space. This 
theory space is spanned by all possible operators which 
can be constructed from the chosen degrees of freedom 
and which are compatible with the desired symmetries. 
If such a fixed point exists the microscopic action to be 
quantized can be constructed from the properties of the 
fixed point and thus is a prediction itself. 

The framework for a nonperturbative construction of 
renormalizable field theories is provided by the functional 
RG which can be formulated in terms of a flow equation 
for the effective average action T)-, the Wetterich equa- 
tion [23[: 

d t r fc [4>] = isTr{[lf ^ + R k ]-i(d t R k )}. (2) 

( 2) 

Here, T k ' [$] is the second functional derivative with re- 
spect to the field the latter representing a collective 
field variable for all bosonic or fermionic degrees of free- 
dom, and Rk denotes a momentum-dependent regulator 
function that suppresses IR modes below a momentum 
scale k. The solution of the Wetterich equation provides 
an RG trajectory in theory space, interpolating between 
the bare action Sa to be quantized Ffc_,A — * S\ and the 
full quantum effective action T = T]~—>o, being the gener- 
ating functional of 1PI correlation functions; for reviews, 
see j24j . 

The effective average action can be parameterized 
by a possibly infinite set of generalized dimensionless cou- 
plings gi. Then, the Wetterich equation provides us with 
the flow of these couplings dtgi = f3 9i (gi, 32, • ■ ■ )• A fixed 
point g* is defined by 

&( 5l *,<7 2 V..) = 0Vi. (3) 
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The fixed point is non-Gaufiian, if at least one fixed-point 
coupling is nonzero g* ^ 0. If the RG trajectory flows 
into a fixed point in the UV, the UV cutoff can safely 
be taken to infinity and the theory can be considered as 
fundamental. 

In addition to being fundamental, we also want the 
theory to be predictive. For this, let us consider the 
fixed-point regime, where the flow can be linearized, 



leading order 



d t g l = Bi 3 (g*-g,j) + 



B 4 > 



d(3 g 



9=9* 



(4) 



The critical exponents Q 1 correspond to the negative of 
the eigenvalues of the stability matrix . They allow 
for a classification of physical parameters: Whereas all 
cigendirections with Q 1 < die out towards the IR and 
thus are irrelevant, all relevant directions with Q 1 > 
increase towards the IR and thus determine the macro- 
scopic physics (for the marginal directions O 7 = 0, it de- 
pends on the higher-order terms in the expansion about 
the fixed point). Hence the number of relevant and 
marginally-relevant directions determines the number of 
physical parameters to be fixed. The theory is predictive 
if this number is finite. In the case of the Gaufiian fixed 
point g* = 0, the present construction corresponds to 
the standard perturbative power-counting analysis and 
the critical exponents are equal to the power counting 
dimensions of the couplings. 

If a critical exponent is much larger than zero, say of 
0(1), the RG trajectory rapidly leaves the fixed-point 
regime towards the IR. Therefore, separating a typical 
UV scale where the system is close to the fixed point from 
the IR scales where, e.g., physical masses are generated 
requires a significant fine-tuning of the initial conditions. 

For the flow towards the IR, the linearized fixed-point 
flow Eq. (|3]) generally is insufficient and the full nonlin- 
ear (3 functions have to be taken into account. Even the 
parameterization of the effective action in terms of the 
same degrees of freedom in the UV and IR might be in- 
appropriate. Nevertheless, we use the same bosonic and 
fermionic degrees of freedom on all scales in the present 
work, since we specifically want to address the question 
whether standard-model IR degrees of freedom can have 
an interacting UV completion. 



III. RENORMALIZATION FLOW OF CHIRAL 
YUKAWA SYSTEMS 

A. Derivative expansion 

In the present work, we investigate a Yukawa theory 
with chiral fermions including one right-handed fermion 
and TVl left-handed fermions. The fermions are coupled 
to TVl complex bosons via a simple Yukawa interaction. 
We span the theory space by a truncated action func- 
tional in a derivative expansion, which reads at next-to- 



d d x{t(Z L , k i'tH't + Z R ^ R # R ) 



(5) 



The fermion fields and ip~R have standard kinetic terms 
but can acquire different wave function renormalizations 
Z^ k and Z-R^k- The index a runs from 1 to N^. The pro- 
jections on the left-/right-handed fermion contributions 
are carried out via the projection operators 



Pl/r=^(1±75). 



(6) 



The bosonic sector involves a standard kinetic term with 
wave function renormalization Z^ k and an effective po- 
tential U k {(j) a ^ <p a ). Defining the invariant p := (f> aji (f> a , the 
effective potential U k (p) can be expanded in powers of p. 
The bosons can also be expressed in terms of a real field 
basis by defining 

</> a = ^W+^2), Qt = ^W-*^), (7) 

where 4>i,4>2 € R. The truncated effective action 
J5|) including the Yukawa interaction is invariant under 
U(-/Vl)l transformations of the left-handed fermion and 
the boson as well as U(1) R transformations of the right- 
handed fermion and the boson. 

All the parameters in the effective average action are 
understood to be scale dependent, which is indicated by 
the index k. 

The flow of the wave function renormalizations 
Z^.k, -^L.fc and Z^ k can be expressed in terms of scale- 
dependent anomalous dimensions 



-dtlnZ^.k, ^3 = -d t lnZ LiRjfc . 



(8) 



Setting the anomalous dimensions to zero defines the 
leading-order derivative expansion. At next-to-lcading 
order, it is important to distinguish between Zl./c and 
Z R k as they acquire different loop contributions, see be- 
low. 

In order to fix the standard RG invariance of field 
rescalings, we define the renormalized fields as 



Z^, ^L,R = ZS^,R. 



(9) 



For the fixed-point search, it is useful to introduce di- 
mensionless renormalized quantities 



P = Z<j>. k k 2 d p, 



Z d>,k Z L,k Z R,k k 



~ 4 hl 



u k {p) = k d Uk(p)\p =k i-2p/Zt, k - 



(10) 

(11) 
(12) 



Detailed information about the derivation of the flow 
equations for this truncation in arbitrary spacetime di- 
mensions d is given in App. [B] For our purposes, we use 
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a linear regulator function R k which is optimized for the 
present truncation [25[ . The flow of the effective poten- 
tial in terms of threshold functions which are given in 
App. [Al reads 

d t u k = -du k + pu' k (d - 2 + 770) (13) 
+ 2v d {(2N L - l)4(u' k ) + l d {u' k + 2pu' k ') 
- d 7 ((iV L - l)zg d (0) + &(~phi) + l$?(phl))}, 

where the primes denote derivatives with respect to p, 
and v d = l/(2 d+1 7r d / 2 r(d/2)). For the symmetric phase 
(SYM), we expand the effective potential around zero 
field, 



Np \ 

E*n,k ~ n 
n! ^ 



A 



m kP+^-p +^rrP +...(14) 



3,fe ~3 



2! 



3! 



For the SSB phase, where the minimum of the effective 
potential u k acquires a nonzero value K k '■= Pram. > 0, 
we use the expansion 



Uk 



\ ~< A„ k . _ 



(15) 



n=2 
^2.k 



! (^-^) 2 +^f(p-- fc ) 3 + 



Given the flow of u k (fill)) , the flows of m 2 or X Uik in 
both phases can be read off from an expansion of the 



flow equation and a comparison of coefficients. For the 
flow of Kfe, we use the fact that the first derivative of u k 
vanishes at the minimum, u' k (K k ) = 0. This implies 



= d t u' k { Kk ) = d t u' k (p)\p =Kk + (d t Kk)u'k(K k ) 

^ d * K h = - -JT7- , dtu' k {p)\~ p=Kk . (16) 

u k \ K k) 



The explicit flow equations for the running parameters 
will be given in the following sections for the SYM and 
the SSB phase separately. Note that the expansion co- 
efficients X n , k in Eqs. (| 14[) and (fTS"]) are not identical. 
Since there is little risk that the notation of the differ- 
ent regimes interferes with each other, we refrain from 
introducing different symbols. 

In the SSB regime, the flow of the Yukawa coupling 
and the scalar anomalous dimension for the Goldstonc 
mode can, in principle, be different from that of the radial 
mode. As the Goldstone modes as such are not present 
in the standard model, we compute the Yukawa coupling 
and the scalar anomalous dimension by projecting the 
flow onto the radial scalar operators in the SSB regime. 
Note that this strategy is different from that used for 
critical phenomena in other Yukawa or bosonic systems, 
where the Goldstonc modes can dominate criticality. 

Accordingly, the flow of the Yukawa coupling h k can 
be derived (see App.[C|), and we end up with 



J 



d t h 



(d - 4 + + r?L + m)h% + 4v d hi{(2pu' k y { ™ )d (phl u' k ) - (6p< + Ap 2 u'l!)lf^ d {phl u' k + 2pu'£) (17) 



-C )d (phl u' k ) + l™(phlu' k + 2-pu'i) + 2phtl^ d (phl u' k ) 2phll^ d (phl u' k + 2pu'l)}. 
Finally, we list the expressions for the anomalous dimensions 

{2N\ J 



Vd, 



^p^ul + 2pu'i'fm d 22 {u' k + 2~pu'l) 



1 ^ Vd :n:" 2 m d r«M 
P u k m 22\ u k) 



8v d d 



\phl) 



^phtrnf^phl), 



8 -^hl[mT d (phl < + 2-pul) + m^ d (phlu' k )], 

^ht[m£ B)d (phl u' k + 2-pu'l) + mT\phlu' k ) + 2(N h l)mg B)d (0, u' k )}. 



(18) 

(19) 
(20) 



r 



The arguments of the threshold functions have to be 
evaluated at the minimum of the effective potential. In 
the following sections, we will concentrate on the case of 
d = 4 dimensions. 

Let us emphasize that the derivative expansion has al- 
ready been tested in various Yukawa systems and has 
proved to yield qualitatively and quantitatively accurate 



results. RG flows for Yukawa systems have been suc- 
cessfully studied in QCD l26l. critical phenomena [27| . 
quantum phase transitions [281 ] and ultra-cold fcrmionic 
atom gases (29j . 
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B. Parameter constraints 



Let us finally discuss several constraints on the cou- 
plings as, e.g., dictated by physical requirements as well 
as by our truncation. As our truncation is based on 
a derivative expansion, satisfactory convergence is ex- 
pected if the higher derivative operators take little in- 
fluence on the flow of the leading-order terms. In the 
present case, the leading-order effective potential receives 
higher-order contributions only through the anomalous 
dimensions. Therefore, convergence of the derivative ex- 
pansion requires 



?/L,??R,^ ^ o{\). 



(21) 



This condition will serve as an important quality criterion 
for our truncation. The SYM regime is characterized by a 
minimum of Uk at vanishing field. A simple consequence 
is that the mass term needs to be positive. Also, the 
potential should be bounded from below, which in the 
polynomial expansion translates into a positive highest 
nonvanishing coefficient, 



"ife, A„ maxjfe > 0. 



(22) 



In the SSB regime, the minimum must be positive, > 
0, the potential should be bounded, and in addition the 
potential at the minimum must have positive curvature, 



K k: ^llmax, k, ^2.k > 0. 



Finally, Osterwalder-Schrader positivity requires 



hi > 0. 



(23) 



(24) 



Beyond that, there are no constraints on the size of the 
couplings as in perturbation theory. 



IV. THE SYMMETRIC REGIME (SYM) 

Let us first investigate the fixed-point structure of the 
system in the symmetric regime. We restrict ourselves to 
four dimensions and expand the effective action using the 
ansatz (| 14[) . Now and in the following we will suppress 
the index k at the parameters of the truncation for no- 
tational convenience. The flow equations are evaluated 
at the minimum of the effective potential, which is at 
p = in the symmetric regime, and we replace u' = m 2 , 
u" = A 2 , u'" = A3, ... . In this case, the next-to-leading 
order flow equations up to second order in the effective 



potential read 

d t m 2 = (770 - 2)m 2 - 



1 



1) 



d t X 2 



2r?0A 2 
1 



1 

8^ 



16tt 2V 6 ' (1 
2 

87T~ 2 V ~ 

6 



(25) 



h 2 . hr VR 



2 (1-^) 



4) 



16tt 2 
An 2 



1 + TO 2 ) 3 



(26) 



(1 



6 



(1 + 77 
h 4 



2\2 



4tt 2 



(1 



5 ' 



We observe that the flow equation of a coupling A„ in- 
volves m 2 , A2, A n+ i. For the Yukawa coupling we ob- 
tain 



d t h 2 = (7^ + t)l + m)h 2 



(27) 



In particular, all contributions from typical vertex trian- 
gle diagrams vanish in the SYM regime. At leading or- 
der in a derivative expansion where r\i = 0, the Yukawa 
coupling does not flow in the SYM regime. At next- 
to-leading order, we also have to take into account the 
anomalous dimensions 



1 

16tt 2 
1 

16tt 2 



1 - 



8 



h 2 (i- v 4 



h 2 (i- v -± 



5/(1 +77Z 2 ) 2 ' 

2Y L 
5 J (1 + m 2 ) 2 



(28) 
(29) 
(30) 



For an interacting fixed point with a non- vanishing h 2 to 
exist, the prcfactor (7/0 + ?/l + ^r) in the flow equation 
for the Yukawa coupling would have to vanish at some 
point in parameter space. This implies that the different 
anomalous dimensions need to have a relative minus sign. 
From Eq. (f2"8")h we read off that a relative sign change 
requires either 77^ > 5 or 77L + 77r > 8. Both options 
would clearly be in conflict with the validity bounds of 
the derivative expansion, demanding at least for 77L r < 
0(1). 

Therefore, a nontrivial fixed point of the Yukawa cou- 
pling in the SYM regime based on the criterion + rjL + 
VR = - even if it existed - would be beyond the reli- 
ability bounds of the derivative expansion. Within the 
validity regime of our truncation, we thus find only the 
Gaufiian fixed point for the Yukawa coupling h* = 0. 
Since the flow of the Yukawa coupling is proportional to 
itself, an initial zero value for h 2 will leave the system 
non-interacting at h? = for all scales. The remaining 
interacting system is purely bosonic and does not show 
any nontrivial fixed point in agreement with triviality of 
4> 4 theory. We conclude that the chiral Yukawa model 
in the SYM regime is not asymptotically safe within the 
part of theory space spanned by our truncation. 
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FIG. 3: Leading-order derivative expansion (local-potential approximation) as well as leading-order polynomial expansion: 
fixed-point values {«* , A|, h* 2 } for the two physically admissible fixed points as a function of TVl. 



V. THE REGIME OF SPONTANEOUS 
SYMMETRY BREAKING (SSB) 

Let us continue our fixed-point search in the broken- 
symmetry regime. The flow equations have a richer struc- 
ture here, since the coupling to the condensate can medi- 
ate further effective interactions. Most importantly, the 
broken regime can support the fixed-point scenario with 
a conformal vev. 



A. Fixed-point search to leading order 

For the fixed-point search, we use a polynomial ex- 
pansion of the effective potential about its minimum, cf. 
Eq. (fT5|) . As a check of the quality of this expansion, the 
convergence properties of physical quantities have to be 
determined with respect to increasing orders in this ex- 



pansion. In general, we expect that this expansion gives 
a good approximation of the full effective potential only 
in the vicinity of the potential minimum. Nevertheless, 
this can still lead to sufficient information for a quantita- 
tive estimate of a variety of physical quantities, as they 
are mostly related to properties of the potential near the 
minimum. 

For a first glance, we use the simplest nontrivial ap- 
proximation of the effective potential, u = ^"{p — k) 2 , 
which facilitates a purely analytical treatment of the 
fixed-point equations. The flow of the effective potential 
(|13p. evaluated around the minimum, then boils down 
to flows for p m in = K and u" — X 2 , with all other cou- 
plings being zero, i.e., v! = 0, and vS n ' = for n > 2. 
We also confine ourselves to the leading-order approxi- 
mation in the derivative expansion, where the anomalous 
dimensions vanish, ?7l,r,0 = 0. The flow equations in this 
approximation read explicitely 



9tK = ~ 2K + 32^ (27Vl - 1} + 32^ (1 + Lx 2 ) 2 ~ h A 2 (l + K h 2 ) 2 ' (31) 

dtX2 = vh {2Nh ~ l)xl + i^(i+ 9 2 A «a 2 )3 ~ i (irW' (32) 



and for the Yukawa coupling 



2 1 h A J / 1 \ 6kX 2 ( 1 2 

dtk ~ lfor* (1 + nh 2 ) \ 2X2K VT+k/? + J ~ (1 + 2kA 2 ) 2 VT+^? + 

1 \ 1/1 1 



1 + nh 2 J (1 + 2kA 2 ) \l + nh 2 1 + 2kA 2 
2nh 2 ( 2 \ 2nh 2 ( 



(l + nh 2 ) \l + nh 2 J (1 + nh 2 ){l + 2n\ 2 ) \l + nh 2 1 + 2kA 2 

I 



The resulting set of fixed-point equations {dfK = 
0, dtX 2 = Q,dth 2 = 0} can be solved analytically. For 
a given N-^. we obtain a large number of fixed points 



but only up to two fixed points fulfill the physical con- 
straints given in section IIII Al For Ni, < 4, we find 
only one admissible non-Gaufiian fixed point (NGFP). 
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FIG. 4: Leading order truncation: Critical exponents as a function of JVl. The fixed point corresponding to the open circles 
has two relevant directions, whereas the fixed point corresponding to the filled circles has only one relevant direction. 



For 4 < Njj < 29, two admissible NGFPs occur. For 
30 < Ni, < 57, there is again only one NGFP, whereas 
larger Ni, do not give rise to any physically admissible 
fixed points. The fixed-point values k*, ft* 2 , and A*, as a 
function of Al can be read off from figure [3] 

The universal critical exponents can be deduced from 
the linearized flow around this fixed point, cf. Eq. (|4|). 

d t g i = B i i(g*-g j )+..., (34) 

where gi = (k, \2,h 2 ) in the present simple truncation. 
The expansion coefficients form the stability matrix, 
the eigenvalues — 9i of which correspond to the critical 
exponents apart from a minus sign. The critical expo- 
nents as a function of Al are shown in Fig. 2] for this 
simple truncation, ft is remarkable that one NGFP has 
two relevant directions and the other even only one rel- 
evant direction with a positive critical exponent. This 
implies that the corresponding model is fixed by only 
two or one physical parameter, respectively. This should 
be compared to the Gaufiian fixed point with perturba- 
tive critical exponents of the order K ~ 2, <d\ 2 ~ 0, and 
0/,. ~ 0, implying that 3 physical parameters need to 
be fixed for Gaufiian models. Moreover, both NGFPs 
exhibit a largest critical exponent being smaller than 
the perturbative maximal value, ma x < 2. As a con- 
sequence, these models defined at the NGFPs have an 
improved hierarchy behavior. 



B. Benchmark fixed point for A^l = 10 

For further investigations and by way of example, we 
concentrate on one particular fixed point. We use a 
fixed point with only one relevant direction, as it is phe- 
nomenologically most appealing due to the reduction of 
physical parameters. We choose Al = 10 in the follow- 
ing, since the corresponding fixed point does not give rise 
to extreme coupling values, implying numerical stability. 
In addition, this particular value of Al may also be of 
relevance for certain unification scenarios. 

Furthermore, we extend the lowest-order polynomial 



expansion studied above to higher operators, cf. Eq. (|15p . 

" p A 

^El^^)". ( 35 ) 

n=2 

where A p labels the higher orders. From a technical view- 
point, the fixed-point analysis becomes much more in- 
volved at higher orders. In particular, we have not found 
any analytical solutions to the higher-order fixed-point 
equations. Due to the intrinsic nonlinearity of the fixed- 
point equations, also a numerical search is nontrivial: a 
complete identification of all possible fixed points appears 
out of reach. 

Therefore, we start with the assumption that the sim- 
ple truncation involving only ft 2 , X 2 and k gives already a 
satisfactory fixed-point estimate. As a next step, we can 
search for a fixed-point of the subsequent higher-order 
system in the vicinity of the previous lower-order fixed 
point. This procedure turns out to be self-consistent and 
can be iterated to higher orders. The results for the as- 
sociated fixed-point values up to A p = 6 are summarized 
in Tab. |TJ The resulting fixed-point values show a sat- 
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N p 


= 2 


55.8 


0.0174 


12.11 












= 4 


56.0 


0.0158 


12.09 


-115 


1.30 ■ 10 4 






N P 


= 6 


57.4 


0.0152 


12.13 


-152 


1.20 ■ 10 4 


-8.76 ■ 10 5 


1.44 ■ 10 8 



TABLE I: Fixed-point values at different orders in the poly- 
nomial expansion of the effective potential. The index A^p 
denotes the largest exponent of p occurring in the polynomial 
expansion (|15p or (|35[). We observe satisfactory convergence 
of the fixed point values. 



isfactory convergence behavior, indicating that the poly- 
nomial expansion is quantitatively reliable already at low 
orders. This convergence property is also confirmed by 
the form of the fixed-point potential in a finite environ- 
ment of the expansion point k, as is depicted in Fig. [5] 
The fixed-point potential of Mexican-hat type remains 
quantitatively stable under the inclusion of higher orders 
- even away from the minimum n. 

Whereas the fixed-point potential is scheme depen- 
dent, also the stability of the universal critical expo- 
nents can be checked under the inclusion of higher or- 
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FIG. 5: Convergence of the fixed-point potential u* upon 
the inclusion of higher orders in the polynomial expansion 
of the effective potential. The dotted line shows the fixed- 
point potential within the simple truncation of Subsect. IV Al 
including only the couplings k and A2. The dashed curve also 
includes A3 and A4, and the solid line extends up to Ae. 



ders in the polynomial expansion. The results for the 
associated critical exponents up to N p = 6 arc summa- 
rized in Tab. [TTJ The leading exponents vary only on 
the 10%-level upon the inclusion of higher-order terms. 
Most importantly, the number of relevant directions is 
not changed at higher orders, and the maximal critical 
exponent tends to smaller values at higher orders which 
leads to an improved hierarchy behavior. We conclude 
that the fixed point exists at leading-order in the deriva- 
tive expansion. It possesses many desirable properties 
and can well be approximated by a polynomial cxpan- 
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FIG. 6: Flow of the leading couplings within the lowest-order 
truncation. 
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= 4 


1.167 
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= 6 


1.056 


-0.175 


-2.35 


-4.97 


-8.49 


-14.02 


-25.54 



TABLE II: Critical exponents in the polynomial expansion 
of the effective potential. The index A?p denotes the largest 
exponent of p occurring in the polynomial expansion (|15[1 or 
(|35[) . We observe a reasonable convergence of the universal 
critical exponents. 



the UV cutoff to 2 /A 2 . 1 

The bare parameters can be fixed in terms of renormal- 
ization conditions. The latter ultimately relates these pa- 
rameters to the physical IR parameters: the Higgs vev v, 
the Higgs mass mniggs and the top mass mt op - These, in 
turn, are related to the renormalized couplings k, h 2 , A2 
by 



v = lim V2k k, 

k^0 



m top = v/1 2 v, mm. 



= v A2 v 



(36) 

(The bottom- type quarks remain massless in our model.) 



VI. PREDICTIVE POWER OF 
ASYMPTOTICALLY-SAFE YUKAWA SYSTEMS 

We illustrate the predictive power of an asymptotically 
safe Higgs sector, using the flow of our system in the 
local-potential approximation as a model in its own right. 

Perturbatively, our model would be defined by three 
bare parameters if the microscopic theory is initiated near 
the Gaufiian fixed point. These three parameters can be 
specified, e.g., in terms of the bare Yukawa coupling h, 
the (f> 4 coupling A2, and the scalar mass term in units of 



This counting of independent physical parameters 
changes drastically if the system originates from the non- 
Gaufiian fixed point in the limit of an infinite UV cutoff. 
As the non-GauBian fixed point has only one relevant di- 



1 Owing to triviality at the GauBian fixed point, the cutoff cannot 
be removed. Strictly speaking, this leads to a lack of universal- 
ity, such that the couplings of higher-dimensional operators also 
correspond to physical parameters. Their influence on the IR 
physics is, however, suppressed by corresponding inverse powers 
of the UV cutoff A. 
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rcction, there is only one physical parameter which we 
relate to the vev v = 246GeV. In fact, fixing v to its 
physical value merely corresponds to fixing the absolute 
scale. Dimensionless couplings or mass ratios are even 
independent of this fixing scale and thus are completely 
parameter-free predictions of the theory. 

In practice, we have to start the flow at some finite 
value of the auxiliary UV cutoff A (in arbitrary units), 
solve the flow, and then associate the resulting value of 
v with its correct value in physical units. This also fixes 
the value of A in units of, say, GeV. In the model with 
Nl = 10, we then end up with a top and a Higgs mass 
at the IR scale 

m top = 5.78d, mHiggs = 0.97t>, (37) 

where the proportionality factors are approached in the 
limit v/A <C 1. In this limit, this result is robust as 
we vary the UV cutoff, which demonstrates the cutoff 
independence of our prediction. As a result, the top and 
the Higgs mass arc pure predictions of the theory itself. 

Whereas the Higgs mass is of the order of the vev as ex- 
pected, mHiggs — 239GeV for v = 246GeV, the top mass 
is significantly larger in our model: m top ~ 1422GcV. 
This is an indirect consequence of the comparatively large 
value of the fixed point of the Yukawa coupling. In or- 
der to elucidate this relation, let us note that the flow 
can be subdivided into three different regimes: first, the 
UV fixed-point regime where all couplings stay close to 
their fixed point values; second, a crossover region where 
all couplings run fast; and third, a freeze-out region in 
the IR where the dimensionful vev builds up and induces 
decoupling and a subsequent approach to the generically 
frozen IR values. 

The dynamics in the crossover region is an intrinsic 
parameter-free property of the model that occurs only in 
a rather narrow window of momentum scales. Hence, the 
couplings have only a finite "RG time" - similar to the 
number of e-foldings in inflationary cosmology - to run 
from their UV initial conditions to their IR- frozen values. 
As the initial fixed-point Yukawa coupling in our model 
is large, a sizeable IR value remains and gives rise to a 
large top mass. The same mechanism has been identified 
in the Z2-symmetric Yukawa system in [ljj even for the 
case with two relevant directions. There, the top mass 
is, in principle, a physical parameter to be fixed; still, 
the finite number of crossover e-foldings leads to a lower 
bound on the Higgs mass. 

Realistic scenarios with a physical value of the top 
mass require either a smaller fixed-point Yukawa cou- 
pling, or a sufficiently large number of crossover e- 
foldings that facilitates a long running of the Yukawa cou- 
pling to its physical value. As the number of e-foldings 
of the crossover window is mainly dictated by the RG 
speed at which the physical vev builds up, large numbers 
of e-foldings can occur if the largest critical exponent 
dominating the k flow is small. We conclude that a re- 
duced hierarchy problem and a realistic value of the top 
mass can go hand in hand with each other. 



Finally, we have to address an artifact of our flow in 
the IR. As can be read off from Fig. [51 the perturbatively 
marginal couplings h 2 and A2 still exhibit a slow log-like 
running in the IR. In our present set of flow equations, 
the source of this running can be traced back to the Gold- 
stone modes which are a particularity of our model, but 
are not present in the standard model. However, even 
if we considered the Goldstone modes as physical, the 
log-like running is still an artifact of our truncation in 
the IR. This artifact arises from the Cartesian decompo- 
sition of the scalar field, <j> = ^5(^1 + £<fe) ( we suppress 
flavor indices in the following argument) where the vev 
is associated with the 4>\ direction and 02 is populated 
with Goldstone modes. This decomposition leads to log- 
like divergencies in the diagrams with internal 4>2 loops 
coupling to external (f>\ legs. However, as becomes clear 
from a proper nonlinear decomposition of the scalar field 
4> = yfpe 16 , the true Goldstone modes encoded in the 
field 6 do not couple to the radial mode via the poten- 
tial, as U{$<t>) = U(p) is independent of 6. We conclude 
that the log-like running would be absent in a proper 
nonlinear field basis in the broken regime [30l ] . 

In practice, we read off our estimates for the IR values 
of h 2 and A2 at a dynamically defined scale km- For 
the identification of this scale, we note that higher-order 
couplings A„>2 are not affected by the Cartesian basis 
and safely flow to zero due to the leading-order Gaufiian 
flow, 

d t X n = 2(n-2)X n + .... (38) 

as expected, see Fig. [7] We define the scale km by that 
scale at which the coupling A3 approaches its zero IR 
fixed point from above within an accuracy of 0.1%. The 
IR predictions listed in Eq. ((37)) have been read off at 
this scale. Of course, in the presence of this flow arti- 
fact, these predictions actually depend on the choice of 
km . Nevertheless, the generic mechanisms in the UV and 
crossover regimes outlined above are not affected by the 
IR artifact. 



VII. PROPERTIES OF THE DERIVATIVE 
EXPANSION AT NEXT-TO-LEADING ORDER 

In the preceeding section, we have analyzed the chiral 
Yukawa model in a derivative expansion of the effective 
action. As a simple criterion for this expansion to be 
valid, the anomalous dimensions measuring the influence 
of higher-derivative terms should be small, cf. Eq. (f2"T]) . 

As a first check of this criterion, we use our fixed-point 
results from the lowest-order polynomial expansion in the 
local-potential approximation and insert them into the 
right-hand sides of the anomalous dimensions Eqs. (|18p ~ 
(f2"0"]) . For this estimate, we ignore the RG improvement 
in the form of the back-reactions of nonzero 77's on the 
flow of the couplings. The results as a function of are 
shown in Fig. [5] Whereas the criterion (|2ip is satisfied 
for rjtj, and 77L for Vl > 5, the anomalous dimension tjr 
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FIG. 7: Flow of the higher-order couplings A3 and A4 of the 
effective potential. The flows of even higher order couplings 
A5 and A5 look similar. 



of the right-handed fermion violates the criterion for all 
N-^. The violation is dramatic for the fixed point with 
two relevant directions and still substantial for the fixed 
point with one relevant direction. 

One reason for the size of tjb, lies in the fact that 
the massless Goldstone modes and massless bottom-type 
fermions contribute strongly. This is because they (i) 
are not damped by massive threshold effects induced by 
couplings to the condensate, and (ii) contribute with a 
large multiplicity ~ N^. We conclude that the leading- 
order derivative expansion does not provide for self- 
consistent estimates of the fixed-point structure of the 
chiral Yukawa model as such. Nevertheless, since the 
Goldstone as well as the massless fermion modes are not 
present in the standard model but a particularity of our 
reduced toy model, we expect that the local-potential 
approximation of the present model provides for a better 
picture of a possible fixed-point structure of the standard- 
model Higgs sector than the literal version of our model 
including its massless modes. 

It is still an interesting question within the chiral 
Yukawa model, whether the nonlinearities induced by 
the anomalous dimension gives rise to admissible fixed 
points if the full RG improvement is taken into account. 
For this, we have first followed the evolution of our 
leading-order fixed points upon gradually switching on 
the anomalous dimensions. Our numerical results are 
compatible with a destabilization of this fixed point, po- 
tentially going along with an annihilation with other fixed 
points. A numerical search of further admissible fixed 
points upon the full inclusion of the anomalous dimen- 
sion remained inconclusive. Owing to the high degree of 
nonlincarity of the coupled equations, a systematic search 
in this high-dimensional parameter space is computation- 



2.5 



2.0 
1.5 
1.0 
0.5 



0.0, 







10 15 20 25 



2.5 
2.0 
1.5 
1.0 
0.5 
0.0, 



ggagaggggjgi 







10 15 20 25 



o°°o. 



100 
50 

20 
10 
5 



°0 



°o, 



00 



°o 



00, 



°o c 



Nl 



5 10 15 20 25 
FIG. 8: Leading-order estimate of the anomalous dimensions. 

ally expensive and beyond the scope of the present work. 

VIII. CONCLUSIONS 

We have explored a new possible route to asymptoti- 
cally safe Yukawa systems with a chiral U(7Vl)l£3>U(1)r, 
symmetry, similar to the Higgs-top sector of the stan- 
dard model. This work is motivated by the fact that an 
asymptotically safe Yukawa sector can solve the triviality 
problem of the standard-model Higgs sector and poten- 
tially improve on the hierarchy problem, both of which 
represent genuine problems related to the UV stability 
and completion of standard-model particle physics. In 
our scenario proposed in [l7| , asymptotic safety is based 
on a conformal Higgs expectation value that induces a 
non-Gaufiian UV fixed point for all essential couplings. 

This conformal behavior can arise from the dynam- 
ics of the theory as a consequence of a balancing be- 
tween fermionic and bosonic fluctuation contributions 
to the running of the vev. More specifically, since 
bosonic fluctuations need to dominate, the existence of 
this conformal-vev behavior depends on the algebraic 
structure of the theory. It is a central result of our work 
that chiral models with an asymmetry between left- and 
right-handed fermions in fact support this boson domi- 
nance. 
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Our quantitative investigations are based on the func- 
tional RG, using the derivative expansion as a nonper- 
turbativc approximation scheme for a systematic calcu- 
lation of the quantum effective action. To leading-order 
(local-potential approximation), our left-right asymmet- 
ric model exhibits the desired non-GauBian fixed points 
for 1 < TVl < 57. Whereas a physically admissible 
non-Gaufiian fixed point establishes the existence of the 
Yukawa model as a fundamental UV complete quantum 
field theory, the physical properties of such a UV com- 
pletion are quantified by the non-negative critical expo- 
nents of this fixed point. Their number corresponds to 
the number of physical parameters, and their value is a 
measure of the naturalness of a large hierarchy (between, 
say, the Planck and the electroweak scale). 

Most importantly, one of the admissible fixed points 
has only one positive critical exponent. This implies that 
only one physical parameter has to be fixed, e.g., the vev 
v = 246GeV, whereas all other IR quantities such as the 
Higgs or the top mass are a pure prediction of the the- 
ory. Moreover, this critical exponent is also smaller than 
the corresponding value in the perturbative domain, such 
that the hierarchy behavior is substantially improved. 

Unfortunately, the fixed point in our toy model is 
destabilized at higher order in the derivative expansion 
due to massless Goldstone and fermion fluctuations. As 
the latter are a particularity of our model and are not 
present in the standard model, we expect the results from 
the local-potential approximation to be of high relevance 
for a more realistic model of the standard-model Higgs 
sector. On the other hand, our results have proved ro- 
bust under an extension of the approximation scheme in 
the boson effective potential. 

As a main result of the present study, a number of 
conclusions about the requirements for a more realis- 
tic model can be drawn. (1.) The appearance of the 
conformal-vev mechanism depends crucially on the alge- 
braic structure of the theory Boson dominance is re- 
quired and can be induced by a chiral left-right asymme- 
try. (2.) If the conformal vev breaks a continuous global 
symmetry, Goldstone bosons are generated. As they are 
not controlled by threshold effects, they may have the 
tendency to destabilize the UV fixed point. A model with 
spontaneous breaking of gauged symmetries will not be 
affected by this problem. (3.) Even though the thresh- 
old effects of the conformal vev are nonperturbative, our 
mechanism does technically not require strong coupling. 
Comparatively small UV fixed-point couplings would also 
be allowed if not phenomenologically preferred. (4.) The 
present model predicts a very large top mass. For a re- 
alistic top mass, either the Yukawa UV fixed-point value 
should not be too large or the crossover from the UV 
regime to the IR freeze-out should extend over many "e- 
foldings", i.e., a larger range of RG time. 

These requirements together with the remarkable pre- 
dictive power of our scenario are a strong motivation to 
extend the present set of ideas to gauged chiral Yukawa 
models, thereby making another step towards a realistic 



standard-model Higgs sector. 
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APPENDIX A: THRESHOLD FUNCTIONS 



The boson cutoff we use is given by: 

yrniy) = (1 - y)0(l - y), 



(Al) 



where y = q 2 /fc 2 , and the fermion regulator rp(y) is cho- 
sen such that y{l + r-o) = j/(l + rp) 2 . Using this regulator 
in the Wetterich equation, we can perform all momen- 
tum integrations analytically. The result can be given in 
terms of threshold functions which read 
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where we have defined r/^ := |(??r + 



APPENDIX B: DERIVATION OF THE 
EFFECTIVE POTENTIAL FLOW 

For the flow of the effective potential, we project the 
Wetterich equation onto constant bosonic fields and van- 
ishing fermionic fields, 



1 



d t U k = ^8Tr{(r% ) +R k )-\d t R k )}\4= cm x.,i > =o, (Bl) 
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where fi is the spacctimc volume. We have to evaluate 
the r.h.s of Eq. (|B1|) . for which we need the T k 2 ^ ma- 
trix. Taking care of the partly Grassmann- valued field 

(2) 

components and the Fourier conventions, T k is derived 
by 



,(2) 
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Adding both parts and using the optimized regulator in- 
troduced in App. [A] we get 



8 t U k = 2v d k c 



-d 7 v d k d 2 
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where the threshold functions arc defined in Eqs. (|A2| ). 
In terms of dimensionless quantities defined in Eq. (|10[) . 
the potential flow turns into Eq. (|13p in the main text. 



^ = (iPl,...,^). 

The transposition refers to flavor as well as Dirac indices. 
For a proper IR rcgularization, a regulator which is di- 
agonal in field space is sufficient and convenient, 



APPENDIX C: DERIVATION OF THE YUKAWA 
COUPLING FLOW 



Rk(q,p) = S(p - q) 



RkB 

R kF 



with a 2Al x 2Al matrix for the bosonic sector 



RkB 



Z^ k 5 ab p z r B 






Z^ k S ab p 2 r B 



cf. App. E]and an (2A L + 2) x (2A L + 2) matrix for the 
fermionic sector 



RkF = - 



z L , k s ab f 






o 










ZrM 






z^r 





\ 



The matrix T k + R k has the same block form as the 
regulator. Therefore, we can evaluate the bosonic and 
fermionic parts separately. We start with the bosonic 
part. Inverting the matrix, multiplying with the deriva- 
tive of the regulator, and taking the supertrace yields 



d t U k B = - 



d d p 
(2n) a 



d t Rk 



2Nt 



1 



Z^^PbIp) 



Z^^PbIp) 



where Pb(p) = P 2 (1+^b(p))- Introducing P F (p) = p 2 (l + 
r F (p)) 2 and d 7 as the dimension of the representation of 
the Dirac algebra, the fermionic contribution reads 



dtUkF 



d d p r d t [Z Ltk r F (p)} 



(2 

(N L -1 



ttY\ 



Z uk {l + r F {p)) 
Z^ k Z^ k P F {p) 



h~ k p + Z L ^ k Z R , k P F (p) 



d t [Zn. k r F (p)} Z h , k Z K . k P F (p) -i 
Z R , k (l + rp(p)) h\p + Z Uk Z R . k P F (p) /' 



For the derivation of the flow of the Yukawa coupling, 
we first separate the bosonic field into a vacuum expec- 
tation value (vev) v and a purely radial deviation from 
the vev, since we are mainly interested in the Yukawa 
coupling between the fermions and the radial mode. 



i 

71 



i 

71 



4>\{p) 
<t>l(p) 



i<t>l(p) 
i4>lip) 



\^(p) 

(v\ 




l 

71 



W 



(p)J 

Va< l (p)7 



setting all A<p2 Goldstone components to zero. The pro- 
jection of the Wetterich equation onto the flow of the 
Yukawa coupling reads 



d t h k 



V25 



2 Sfoip) 6A</>{(p>) 4 k 'sMq) 



(CI) 



The vertical line indicates that the equation is evaluated 
at = ipl = A0 = 0, p' = p = q = 0. Next, we 

(2) 

can decompose the matrix (F^. + R k ) into two parts. 

(2) 

One part, which we call (T k ' Q + R k ), contains only v and 
is independent of the fluctuations. The remaining part, 

(2) 

Ari , contains all fluctuating fields. Inserting this into 



equation (|C1[) and expanding the logarithm 
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STr 



(ln(4 2) +J R fe )) = STr 



In 



, (2 . „ a Ar 



(2) 



Rk 



(C2) 



STr (ln(r^ + J? fc ) +STr 



(2) 



r 



(2) 
fc.O 



Rk 



-STr 



Ar 



(2) 



-STr 



Ar 



(2) 



only the term to third power in Arjf 2 ' survives the projection onto A^JV'lV'R- Performing the matrix calculations 
and taking the supertrace, we get 



d t hk 



d d P - 
(2ir) d * 



Z htk Z Rik P F (p) 



3U',!v + U'l'v 3 



(Zt, k p B (p) + u' k y (z^Pnip) + u' k + u^y 



(Z L , k Zn,kP F (p) + f v 2 ) 2 \Z*,kPa{p) + U' h Z^ k P B (p) + U' k + U£v 
1/1 1 



where the potential on the right hand side is evaluated at the minimum \v 2 . Using the optimized regulator and 
the threshold functions as defined in App. \^ and switching over to dimcnsionlcss quantities, we end up with the 
representation (|17|) given in the main text. 

I 



APPENDIX D: DERIVATION OF THE 
ANOMALOUS DIMENSIONS 

Por the derivation of the flow of Z^^, we decompose 
the bosonic field as in App. [C] The projection of the 
Wetterich equation onto the boson kinetic term leads us 
to 



Again, we have decomposed the matrix into two parts, 
one part containing the vev and the other part contain- 
ing all fluctuating fields. As done in App. [Cl we have 
expanded the logarithm as in Eq. (|C3j) . but this time 
only the second order contributes. Calculating the r.h.s 
of Eq. (|D1[) . results in 



d t Z, 



8 



d{p' 2 ) 6A(f>{(p>) 5A<j)\(q>) 



(Dl) 



x — STr 

4 



Ar 



.2) 



,(2) 



R k 



A(j>=lf>l=lpB_=0,p'=g'=0 



d t Z, 



i f d d P 5 r 

dj (2nr 0t 



+ (2N L - \){UlvYp*Z^ k 



(3UZv + U£'v*)VzIh 



dp 



9 Pb( P ) 



{z^ k Pv{p) + u' k + u' k ^Y 

\ 2 



+2hld^p 4 Z htk Z^ k 







{Z^ k P B {p) + U' k f 
(l+r F (p)) 



d P 2 Z, 



- hld^v 2 p 2 



L,kZ R , k P F (p) + f v 2 J \ d P 2 Z L , k Z R , k P F (p) + 



where the potential is evaluated at the minimum. Using the optimized regulator together with the threshold functions 
as defined in App. [A] and again introducing dimensionless quantities, we end up with 



p(3u'l + 2~pu' k ") 2 m d 22 {u' k + 2~pu' k r ) + 



{2N L - l)8v d x _ „ 2 



pufmi 2 (u' k ) 



-^hlmT\ph 2 k ) - ^phWP\ph 2 k ), 
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where we have used r]^ 



and the right hand side has to be evaluated on the vev. 



For the fcrmionic anomalous dimensions, the procedure is the same. We start with 



1 u d 
tr7 M 



STr 



Ar 



(2) 



, L k,0 +K k, 



A(p=lp- L =tpB.=0,p' =q' =0 



(D2) 



and get 



2hl f d d p 9 ~ 

1 



Z L , k (l + r F (p)) d 



(z^ k p B ( P ) + u' k + u^y 



(Z^ k P B ( P ) + U' k r 



and 



dtZ 



't^R.fc 



d d P 2 , 



(27T) 
1 



z r mi + Mp)) d 

, 2(iV L -l)Z Rifc (l + r F (j5)) 



(9 



(^, fc P B (p) + ^ + ^ 2 ) 2 
1 



{Z^ k P B {p) + U' k f ' Z L , fc ZR, fc iV(p) "^V^^fePBW + C/D 2 



respectively. Using the optimized regulator the result read in terms of dimensionlcss quantities 



>1L 



hl[m^ d (phl,u' k + 2p<) + ^"^.ui)], 



(D3) 



and 



m = ^hi[m£ B)d (phl,u> k + 2p<) + m^ d (phl,u' k ) + 2(N h - l)m^ B)d (0,<)]. 



(D4) 
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